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ABSTRACT 

We revisit the limits of the level of the extragalactic background light (EBL) recently reported 
by the MAGIC collaboration based on the observed 7-ray spectrum of the quasar 3C279, con- 
sidering the impact of absorption of high-energy 7-ray photons inside the broad line region 
(BLR) of the quasar We use the photoionization code CLOUDY to calculate the expected 
optical-UV radiation field inside the BLR and the optical depth to 7-rays for a relatively ex- 
tended set of the parameters. We found that the absorption of 7-ray photons, though important 
for the estimate of the true radiative output of the source, does not produce an important hard- 
ening of the spectrum of 3C279 in the energy band accessible by MAGIC, supporting the 
method used to infer the upper limits to the level of the EBL. 
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1 INTRODUCTION 

The detection of 3C279 (2;=0.536) in the very high energy (VHE, 
S > 50 GeV) band by the MAGIC telescope (Albert et al. 2008) 
extends to the quasars the group of the known extragalactic VHE 
sources, before limited to BL Lac objects (excluding the nearby 
radiogalaxy M8713. 

For some aspects, the detection of 3C279 comes as a surprise. 
General theoretical arguments support the view that quasars cannot 
be important VHE emitters, in particular because of the expected 
absorption, through the pair production process 7 + 7 — > + e~ , 
inside the source itself (among the most recent calculations, e.g., 
Donea & Protheroe 2003, Liu & Bai 2006, Reimer 2007). More- 
over, quasars are generally located at a relatively high redshift, 
implying a huge absorption by the extragalactic background light 
(EBL, e.g. Kneiske et al. 2004, Primack et al. 2005, Stecker et al. 
2006, Franceschini et al. 2008). Last, but not least, widely adopted 
standard leptonic models for production of 7-rays in 3C279 (Hart- 
man et al. 2001, Ballo et al. 2002; see also Tavecchio & Ghisellini 
2008, hereafter TG08) do not predict an important emission above 
few tens of GeV, because of the rapid decrease of the scattering 
cross section (for hadronic scenarios see, e.g. Mannheim 1993). 
The observation of VHE photons from 3C279 by MAGIC (Albert 
et al. 2008) demonstrates that also quasars can, to some extent, pro- 
duce high-energy 7-rays and suggests that the opacity (both intrin- 
sic and cosmic) is less strong than what previously assumed. 

The detection of a source of VHE photons at a relatively high 



redshift offers a unique tool to probe the still poorly known EBL. 
Albert et al. (2008) used general arguments based on a limiting 
spectral slope for the emitted spectrum to infer the amount of ex- 
tragalactic absorptior0. However, as recently pointed out by Aha- 
ronian et al. (2008) (see also Bednarek 1997), intrinsic absorption 
of 7-ray photons inside the source could result in rather hard ob- 
served spectra. Such spectra affected by absorption could lead to 
severely underestimate the level of the EBL when the arguments 
based on the hardness of the spectrum are used. For this reason, the 
use of the measured spectrum to constrain the EBL in Albert et al. 
(2008) triggered some discussion on the role and strength of the in- 
ternal absorption in this source (Liu et al. 2008, Sitarek & Bednarek 
2008). 

We emphasize that the absorbed spectrum will be harder than 
the intrinsic one only in the case of an optical depth t{E) decreas- 
ing with energy, E. Thus, as long as the optical depth increases (or, 
at least, only slightly decreases) with energy, the resulting spectrum 
will be softer (or slightly harder) than the intrinsic one and the stan- 
dard spectral methods to constrain the EBL could be safely used. 
It is thus clear that the shape of t{E), strictly related to the spec- 
trum of the target photons, is the key element to assess the effects 
of absorption on these methods. 

Previous attempts to calculate intrinsic absorption in the BLR 
(Liu & Bai 2006, Reimer 2007, Liu et al. 2008) assumed rather 
idealized templates for the BLR spectrum, considering the most 
prominent emission lines but neglecting the important contribution 
of the optical-UV continuum. In a recent paper, Sitarek & Bednarek 
(2008) include the absorption in a self-consistent model for the 
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Figure 1. Sketch of the geometry assumed in the model (not to scale). The 
(uniform) BLR is assumed to be a spherical shell with thickness AR and 
inner radius Rblr, illuminated by the central continuum with luminosity 
Lq. The source of 7-ray photons is at distance x from the central black 
hole. See text for details. 

high-energy emission of 3C279. They consider in detail the geom- 
etry of the radiation fields inside the BLR, including also the con- 
tribution of direct and scattered radiation of the accretion disc, but, 
again, they model the BLR radiation with the simplified template 
of Liu & Bai (2006). 

In this work (Sect. 2) we explore the effects of absorption 
using more realistic spectra of the BLR calculated with the pho- 
toionization code CLOUDY (Ferland et al. 1998), previously used to 
study in detail the inverse Compton emission from powerful blazars 
(TG08). In particular, we show that the IR-optical-UV continuum 
plays an important role in determining the absorption, resulting in 
an almost constant optical depth in the broad energy range 30 GeV- 
30 TeV. In Sect. 3 we use the spectra corrected for absorption to 
revisit the constraints on the EBL based on the spectrum of 3C279 
derived in Albert et al. (2008). In Sect. 4 we discuss the results. 



2 INTRINSIC ABSORPTION 
2.1 The model 

We calculate the diffuse radiation field inside the BLR of 3C279 
for different values of the BLR radius, temperature of the accretion 
disk, slope of the illuminating UV radiation. We refer to TG08 for 
a full description of the model. We assume the geometry shown in 
Fig[T] The accretion flow illuminates the BLR clouds (character- 
ized by the total hydrogen density n and the hydrogen column den- 
sity Nh) isotropically filling the BLR, assumed to be a thin spheri- 
cal shell with inner radius -Rblr- In the calculation we assume that 
the clouds cover a fraction C = f2/47r = 0.1 of the solid angle 
viewed from the central illuminating source. The emission from 
the illuminated face of the clouds is calculated with version 05.07 
of CLOUDY, described by Ferland et al. (19980 For simplicity, we 
discuss only the case of solar abundance and in all calculations we 
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fix n = 10"' cm~^ and Nh = 10^^ cm"^. Results do not sub- 
stantially change for different densities and column densities (see 
TG08). We adopt the spectrum of the illuminating continuum mod- 
eled as a combination of a UV bump (with slope auv) with a flat 
X-ray power-law, Ld{v) (x u^^, commonly assumed in these cal- 
culations (AGN model in CLOUDY, e.g. Korista & Goad 2001 and 
references therein). When not explicitly noted, the disk is assumed 
to have the "standard" temperature Td = 1.5 x 10'' K. 

We assume that the disk in 3C279 emits a total luminos- 
ity Ld = 2 X 10*^ erg/s (Pian et al. 1999). More uncertain is 
the radius of the BLR. The empirical relations connecting the lu- 
minosity of the disk and -Rblr (Bentz et al. 2006, Kaspi et al. 

2007) provide -Rblr — 1 — 3 x 10^^ cm. Given these uncertain- 
ties, below we show the results for four different values of -Rblr, 
1, 1.6, 3.2 and 6.3 x lO" cm. 

The optical depth for the photon-photon absorption, t{E), is 
calculated as (e.g., Liu & Bai 2006): 

where E is the energy of the 7-rays, I is the distance of the pho- 
tons from the BH, n(i', Q, I) is the number density of the radiation 
for solid angle at each location of the photon path, a-y-y {E, v, Q.) is 
the cross section and dSl — 2ndfi. Note that, given the characteris- 
tics of the cross section, for a fixed geometry t{E) oc uniy) with 
v oc 1/E, that is the optical depth reflects the spectrum of the soft 
photon field (smeared by the cross section). 

To calculate t[E) from Eq.Q one has to specify the loca- 
tion of the source, x. All calculations discussed below have been 
performed with s = 2 x lO^'' cm. Larger values of x (a source 
closer to the BLR) result in a lower level of the absorption, but 
do not essentially affect the shape of t{E). For smaller values of 
X one should also consider the absorption induced by photons di- 
rectly coming from the accretion disk (Ghisellini & Madau 1996, 
Sitarek & Bednarek 2008). However, absorption from the disk is 
characterized, for the considered range of energies, by an optical 
depth monotonically increasing with energy (Sitarek & Bednarek 

2008) . As already discussed, in this conditions the resulting spec- 
trum is softer than the emitted spectrum and thus it does not affect 
the constraints on the EBL. 

Note that, for simplicity, we neglect the effects related to the 
radiative transfer inside the emission region (negligible as long as 
the size of the source is much less than -Rblr — a::)- Another pos- 
sible effect that we do not take into account is that in the case of a 
moving source (as in the standard "internal shock" scenario, Spada 
et al. 2001), photons emitted at different times will also be charac- 
terized by different x in Eq[T]and thus will suffer a different level of 
intrinsic absorption. The study of these effects, in part already con- 
sidered by Sitarek & Bednarek (2008), is important in view of the 
interpretation of the time-resolved spectra soon available thanks to 
the Fermi Gamma-ray Space Telescope. 



2.2 Results 

As an example, Fig|2| reports some BLR spectra derived with the 
model, assuming Td ~ 1.5 x 10^ K, auv ~ 0.5, -Rblr, 
1, 1.6, 3.2 and 6.3 x 10^^ cm, plotted as the number photon 
density, n{i')i'. The two vertical lines indicate the spectral range 
interesting for absorption of 7-rays in the energy range covered 
by the MAGIC spectrum (upper x-axis, in the quasar rest frame). 
Clearly, besides the emission lines, the continuum (deriving from 
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Figure 2. Number density of BLR photons for = 2 X 10^^ erg/s, 
Td = 1.5 X 10^ K, auv = 0.5 and BLR radii of I, 1.6, 3.2 and 
6.3 X 10^^ cm (cases S in Fig. 3). Note tlie important component in the 
optical-UV band, mainly coming from from free-free emission, Thomson 
and Rayleigh scattering. The upper x— axis reports the energy of 7-rays 
mainly absorbed by soft photons of the corresponding frequency. The two 
dotted vertical lines show the frequency range (in the quasar rest frame) in- 
teresting for absorption of 7-rays in the interval of energies covered by the 
MAGIC spectrum reported by Albert et al. (2008). 



free-free emission, recombinations and scattering, both Thomson 
and Reyleigh) provides an important contribution especially at low 
frequencies, below the Lya line. In particular, in the UV band the 
dominant contribution comes from Thomson and Rayleigh scatter- 
ing (e.g. Korista & Ferland 1998), whereas at IR-optical frequen- 
cies (in the range v = 10^^ - 10^* Hz), the emission is generally 
dominated by optically thin free-free emission plus a contribution 
from Thomson scattering. The result is that the n{v)u curves are 
almost flat below the UV band, translating in an almost constant, 
or slightly increasing, optical depth for a rather broad interval of 
energies. 

Fig[3] illustrates some examples of the optical depth, t{E) 
(upper panel) and the corresponding absorbed spectra, assum- 
ing an intrinsic photon spectrum F{E) oc _E~'^'^ (lower panel). 
Different line styles and colours refer to different parameters. 
Short dashed, long dashed, dotted and solid lines are for 7?blr, 
1, 1.6,3.2 and 6.3 x lO^'^ cm, respectively. Black lines refer to 
the "standard" (S) scenario, with auv ~ 0.5 (Elvis et al. 1994) 
and Td = 1.5 x 10^ K. Red lines ("extreme" case, E) are calcu- 
lated for a somewhat extreme value of the UV slope, auv = — 1/3 
(slope expected from a standard thin disk), while green lines report 
the results for a "low temperature" case (L), Td = 5 x 10* K (and 
auv = 0.5). 

A common feature of these curves is the sudden increase of 
t{E) starting around i5 ~ 20 — 30 GeV, energy at the threshold for 
photons of the Lya line. After reaching the maximum around 100 
GeV, ''"(£') displays a decreasing branch, just centered on the en- 



Log V [Hz] 
16 15 14 13 12 




1 2 3 4 5 6 
Log E [GeV] 



Figure 3. Upper panel: optical depth for absorption of 7-rays for different 
parameters. Short dashed, long dashed, dotted and solid lines: -Rblr = 
1, 1.6, 3.2 and 6.3 X lO^'^ cm, respectively. Black lines (case S): O-uv = 
0.5andrD = 1.5x10^ K.Redlines (caseE):a(7y = -1/3. Green lines 
(case L): Td = 5 x 10* K. The two vertical lines indicate the energy range 
of the observed spectrum of 3C279 (quasar rest frame). The upper x-axis 
shows the frequency of the target photons mainly interacting with 7-rays of 
the energy indicated in the lower a;-axis. Lower panel: modification of the 
intrinsic photon spectrum (assumed to be a power-law with slope F = 1.5, 
black dashed-dotted line) by absoiption. For simplicity we do not report the 
curves of the L case, always leading to a softening of the spectrum in the 
MAGIC band. Lines styles and colours as above. 

ergy range covered by the MAGIC spectrum (vertical lines). With- 
out the contribution of the continuum at frequencies below the Lya 
line, the optical depth would fast decrease, determining a very hard 
observed spectrum above 100 GeV (rest frame), as in the model of 
Aharonian et al. (2008). However, the important contribution of the 
continuum produces a bump extending above the TeV band, ham- 
pering an important hardening of the out-coming 7-ray spectrum 
(FigEl lower panel). 

Some general characteristics of the optical depth in the "line- 
dominated" (below 100 GeV) and the "continuum dominated" 
(above 100 GeV) regions are clearly evident in these curves. The 
"line-dominated" bump is strongly depressed in the L case. This 
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Figure 4. Energy spectrum of 3C279. The measured spectrum (open circles), EBL-corrected (filled black circles) and four different scenarios for the internal 
absorption (in addition to the external absorption) are shown. 



is due to the small fraction of photons above the Lya energy, with 
the subsequent depression of the luminosity of the Lya line. In this 
case the optical depth is always increasing with energy in the in- 
teresting energy range. Another trend is visible by comparing the 
curves corresponding to the S and the E case: in the latter the 'line- 
dominated" component is systematically more prominent than in 
the S case. The reason is that, for the same luminosity, the hard il- 
luminating continuum in the E case has a larger fraction of photons 
above the Lya energy than in the S case. A third trend visible in 
Fig|3] especially for the L and the S curves, is the increasing role of 
the "continuum dominated" region when J?blr decreases. The rea- 
son of this eiTect is the increasing number of free electrons available 
for scattering and free-free emission due to the increasing ioniza- 
tion at small radii. For the E case the importance of this effects is 
partly reduced, because of the paucity of the small optical-IR flux 
of the illuminator. 

As we have already stressed, the hardening of the spectrum 
is realized for decreasing values of the optical depth with energy. 
The key parameter determining the slope of the optical depth (and 
thus the slope of the modified 7-ray spectrum) at energies above 
~ 100 GeV is the ratio between the optical depth in the "line- 
dominated" bump and in the "continuum dominated" region. This 
is a robust consequence of the realistic models of the BLR radiation 
considered here. Therefore the hardest spectra will be observed in 
the cases of the largest value of the ratio line/continuum. As dis- 
cussed above these conditions are realized for large temperatures 
of the disk and hard UV slopes. Case E can then be considered a 
conservative upper limit for the calculations. 

Besides the spectral modifications, an important aspect to con- 
sider is also the level of the absorption. In general, smaller radii 
imply larger absorption, since r oc 1 / R. In the case of the smaller 
radius considered here, Rblr = 10^^ cm, fluxes are depressed by 
more than 2 orders of magnitude for the S and the E cases, pushing 
the power requirements of the source to above 10^" erg/s. The ab- 



sorption in the MAGIC band is instead limited in the case of the L 
case. For the largest radius, Rblr ~ 6.3 x 10^^ cm, the require- 
ments are still large for the E case (r ~ 2), while for the other two 
cases the absorption is modest. 



3 LIMITS TO THE EBL 

In the following we discuss the effect of absorption on the con- 
straints for the EBL. We use the EBL model from Kneiske et al. 
(2004), modified in Albert et al. (2008) to represent an upper limit 
of the EBL level, which is in the same time a lower limit on the 
transparency of the universe to VHE 7-rays. Using this particular 
EBL model we examine different scenarios of the internal absorp- 
tion discussed above and revisit a possibility to emit corresponding 
VHE spectra. From the discussion above it is clear that the modi- 
fication of the spectrum in the E case should be considered a con- 
servative upper limit to the real case. The L case is not considered 
since it always leads to softer spectra. Moreover we conservatively 
consider only the cases for which the ratio of the absorbed and the 
emitted flux in the MAGIC band is larger than 10~^ (corresponding 
to r < 4.5). 

For the discussed scenarios of the internal absorption, we re- 
constructed the intrinsic spectrum of 3C279. First, we corrected 
the measured energy spectrum (Albert et al. 2008) for a given EBL 
model, corresponding to the maximum allowed level derived in that 
paper. In such a way we reconstruct the energy spectrum which es- 
capes the vicinity of 3C279 towards the observer In the next step, 
we corrected this spectrum for the intrinsic absorption to obtain the 
original (produced) energy spectrum. Results are shown in Fig|4] 
As discussed above, it can be seen that different scenarios for the 
internal absorption mainly affect the flux level of the emission but 
not the shape of the spectrum. Adopting the test prescription from 
Mazin & Raue (2007), we tested the resulting intrinsic spectra for 
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the criterion of the hardness (F > 1.5) and found that in all tested 
cases the intrinsic spectra can be excluded. Consequently, identical 
or even harder EBL limits can be derived as compared to the ones 
obtained in Albert et al. (2008). Note that for the tests not only the 
fit value of the slope from a simple power law (PL) was examined 
but also slope values from a fit by a broken power law (BPL). The 
latter one was tested in case the BPL fit gave a significantly smaller 
valuflthanthePL. 



4 DISCUSSION 

The detection of VHE 7-rays from 3C279 and EBL constraints de- 
rived in Albert et al. (2008) triggered some discussion on the role 
and strength of the internal absorption in this object (Liu et al. 2008, 
Sitarek & Bednarek 2008). In this paper, we have calculated vari- 
ous models for the BLR diffuse radiation and applied them for the 
case of 3C279. We used the code CLOUDY to calculate the mod- 
els in order to have a realistic energy spectrum of the photons in- 
side of the BLR region. The main difference between our model 
and those of Liu et al. (2008) and Sitarek & Bednarek (2008) is 
the assumed BLR spectrum. Liu et al. (2008) assume that the BLR 
spectrum consists only of narrow lines. Sitarek & Bednarek (2008) 
present a detailed model for the high-energy emission of 3C279, 
including also the possible contribution of direct and scattered disc 
radiation. However, their spectrum of the BLR emission does not 
include other important contributions to the continuum, such as the 
free-free emission. As we have shown, our model provides a strong 
continuum component, extending on the optical-UV band (TG08). 

Our results confirm that radiation from the BLR modifies the 
primary emission of 3C279. However, we find that the internal ab- 
sorption inside the BLR does not produce an important hardening 
of the spectrum in the energy band covered by the MAGIC observa- 
tion. In particular, we found that for the tested BLR models, despite 
a possible overall softening of the 3C279 spectrum, at least part of 
the spectrum was significantly above an implied maximum hard- 
ness of r — 1.5, confirming the EBL constraints derived in Albert 
et al. (2008). Of course, the consideration of intrinsic absorption 
implies that the TeV emission from 3C279 could be substantially 
more powerful than published. 

Note also that conclusions from Sitarek & Bednarek (2008) 
that an EBL model of Stecker et al. (2006) does not imply an un- 
realistic intrinsic spectra of 3C279 (when taking into account inter- 
nal absorption in the BLR) concern the "baseline" EBL model of 
Stecker et al. (2006). The EBL limits in this paper and in Albert et 
al. (2008) instead concern the "fast evolution" model of Stecker et 
al. (2006), implying a significantly higher EBL level in the redshift 
range between z = and z — 1. 

We finally note that the discussion on the role of absorption 
implicitly assumes that the highly variable high-energy 7-ray emis- 
sion detected from 3C279 is produced internally to the BLR (prob- 
ably through the comptonization of the BLR photons). However, 
given the small size of the BLR in 3C279 as estimated from the 
empirical relations connecting it to the disk luminosity, it is con- 
ceivable that the emission is (at least in some occasions) produced 
outside the BLR, thus avoiding the problems connected to absorp- 
tion. In this case, the mechanism responsible for the production 
of the observed emission cannot be the external Compton: alter- 
natives include synchrotron self-Compton emission (Lindfors et al. 

^ criterion was defined using the likelihood ratio test, for details see Mazin 
& Raue (2007) 



2006) or comptonization of the IR radiation from the putative dusty 
torus surrounding the central regions (e.g. Sikora et al. 2002, 2008; 
Sokolov & Marscher 2005). 
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